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Abstract. - We study electron molecules in realistic vertically coupled quantum dots in a 
strong magnetic field. Computing the energy spectrum, pair correlation functions, and dy- 
namical form factor as a function of inter-dot coupling via diagonalization of the many-body 
Hamiltonian, we identify structural transitions between different phases, some of which do not 
have a classical counterpart. The calculated Raman cross section shows how such phases can 
be experimentally singled out. 



Electron systems form a Wigner crystal at sufficiently low density or high magnetic field 
B [0. Theoretical [|||| and experimental studies suggest that lowering dimensionality 
favors localization: in this perspective interacting electrons confined in a quantum dot (QD) 
||, sometimes called Wigner molecules 0], are interesting in their own right J|,|§], due to the 
interplay between the electron-electron repulsion and the confining potential. This leads to a 
complex zero temperature phase diagram [[l0| , as compared to the infinite layer well 
as to complex melting mechanisms M. The formation of coupled QDs (artificial molecules) 
introduces qualitatively new physics jfOJ. New energy scales appear — inter-dot tunneling, 
inter- vs intra-dot Coulomb correlation — , whose balance controls the phase diagram ]l2|] . 
Significantly, these parameters can be tuned by inter-dot distance d and/or electron density, so 
that the nature of these few-particle systems and their phases can be explored experimentally. 

In this Letter we discuss quantum mechanical calculations of N electrons in a coupled QD 
structure in the strong field regime, where localization is ensured in the parent isolated QDs. 
Monitoring the spatial correlation functions, we identify different ground states depending 
on the inter-dot coupling. More interestingly, we find that some of the phases do not have 
a classical counterpart p3| , and are ascribed to a two-three-two dimensional (2D-3D-2D) 
transition of the electronic system. Such phases were not identified in previous studies of the 
coupled layer system due to the neglect of tunneling and/or finite width of the layers ]l4|-|l7|], 
i.e., of the 3D nature of the system, which turns out to be essential for their formation. From 
the analysis of the dynamical form factor, we associate to each phase peculiar collective charge 
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excitations. By explicit calculation of the optical spectra, we predict that Raman selection 
rules and peak positions may clearly discriminate between the different phases. 

We consider two vertically coupled QDs, defined by the potential V(r) = to*(jJq£ 2 /2 + 
V(z) given by the sum of an in-plane term [g = (x,y), m* effective mass, loq characteristic 
frequency], and a profile V(z) along the growth direction z, which is a symmetrical square 
double quantum well (DQW). Each well, of width L w and potential height Vo, contains one of 
the two QDs. They are separated by a barrier of width d. The single-particle Hamiltonian in 
the symmetric gauge is Hq(t, s z ) — (—iKV + \e\ A/c) 2 /2m* + V(r) + g*/iBBs z , with B parallel 
to the z-axis (A = B x g/2, hb Bohr magneton, g* effective giromagnetic factor, s z = ±1/2 
spin). Since we are interested in the strong localization regime of small filling factors is, 
we expect that spin texture does not alter the essential physics. Therefore we assume that 
electrons are spin polarized Es| and neglect Landau level mixing The eigenfunctions 

of Hq are ij) m ,g( r ) = tPmio) Xi( z )i where ip m {_Q) (m = 0, 1, 2, . . .) are the Fock-Darwin orbitals 
of the first Landau band, and Xi and Xi the symmetrical (S, bonding) and antisymmetrical 
(AS, antibonding) states of the DQW, respectively. We neglect higher subbands since in real 
QDs the confinement in the z direction is stricter than in the xy plane. The index g in ip m ,g(r) 
stands for the parity under spatial inversion r — > — r: even when m + i is odd, odd otherwise. 
The many-body Hamiltonian H. [a = (to, g)] is: 

H = J2 £ o c Ua + 1^2 (H K r |ri-r 2 | l 7<5 ) c ° c ^ C7Ci - 

c a destroys an electron occupying the orbital a. Here the single particle energy e mg — hCl(m + 
1) — huj c m/2 + 6i — \g*\fiBB /2 is the sum of the in-plane contribution and the energy of the 
z-th DQW state e,, f2 = (ljq + w^/4) 1 / 2 , uj c — \e\B/cm* is the cyclotron frequency, n r is the 
dielectric constant |l- 7i is represented in a basis of Slater determinants spanned by filling 
with N electrons the single-particle states ip m ,g] h is diagonalized on each Hilbert space sector 
labeled by the total orbital angular momentum M and parity. 

Since the basis must be truncated, as in any Configuration Interaction approach, the Fock 
space was built by filling up to 32 orbitals, chosen within the set {tp m ,g}m,g to minimize the 
energy. In order to improve the accuracy of results, the choice of single-particle orbitals ipm,g 
depended on the M-sector: by trial-and-error, we found two optimized sets of orbitals: for 
M > 70, we included S levels with m = 0, . . . , 24 and AS levels with m = 0, . . . , 6, while for 
M < 69 we included S levels with m = 3, . . . , 17 and AS levels with rn = 0, . . . , 16. Subspaces 
obtained in this way (with maximum size « 2 ■ 10 4 ) were diagonalized via Lanczos method. 
The two-body Coulomb matrix elements of TL were computed numerically. 

We consider B fields such that single-QD correlation functions show strong localization J]l| . 
We present results for N — 6. The system with N < 6 exhibits the same physics. We can 
identify three regimes which, in general, correspond to different electron arrangements: (I) 
At small d, tunneling dominates and the system behaves as a unique coherent system. (II) 
As d is increased, all energy scales become comparable. (Ill) When eventually tunneling is 
suppressed, only the ratio between intra- and inter-dot interaction is the relevant parameter 
for the now well separated QDs. 

Figure |l| shows the calculated ground-state energy vs d. At small distances the curve 
increases with d (phase I), because the kinetic energy exponentially grows due to the progres- 
sive localization of the wavefunction into the dots: electrons occupy only S orbitals, whose 
energies increase. Close to d = 5 nm (phase II), inter-dot tunneling is suppressed, so that S 
and AS orbitals are almost degenerate: the dominant energy contribution in this regime is 
the inter-dot Coulomb interaction oc 1/d (phase III). As d varies, ground- and first-excited 
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states, labelled by different M's and parities, cross each other at several critical values (thin 
vertical lines). The monotonous decrement of M jlij proceeds by steps of 5 and then 3 units, 
in regions I and III, respectively; as discussed later, these are "magic numbers" 0, reflected 
in the spectrum of collective excitations. Correspondingly, the in-plane average radius (g) 
(right axis in Fig. has a stair-like behavior and it is almost constant at a given M. M 
measures the in-plane spatial extension of the charge density: the higher M the outer orbitals 
occupied [^0|. We define a total filling factor v, in analogy with double layers in the fractional 
quantum Hall effect (FQHE), as v = N(N - 1)/2M @. Here v ranges from 1/5 at d = to 
5/12 when the two dots are isolated (then v/2 — 5/24 refers to a single dot with N = 3). 

The flat steps of (g) in Fig. |l] imply that these states are incompressible, in the same 
sense as Laughlin's states of the FQHE plfl . Indeed, varying d acts like an external pressure 
applied in the z direction, forcing the wavefunction to change: however, due to a cusp-like 
structure of the energy spectrum |pfj| , this happens only in a discontinuous way. In the 
tunneling-dominated regime, up to d — 4.5 nm (phase I), increasing d means enlarging the 
volume of basically a unique QD, thus forcing a rearrangement of the incompressible charge 
density at the M-transition 75 — > 70. In the "compressible" region II the penetration of the 
wavefunction into the inter-dot barrier goes to zero: the effective volume of the charge density 
is thus reduced, which explains the slight and continuous increase of (g) with d, as shown in 
Fig. 0. For d > 5.3 nm (phase III), the dots are well separated, and v increases from 1/4 to 
5/12 when d — > oo, which is well known from double-layer physics |||2^], due to decrease of 
inter-dot correlation (stabilizing the Wigner molecule). 

To analyze the internal structure of the molecule, we compute the pair correlation function 
P(g, z; q , z ) = ( s (e ~ Qi)S{z - z l )8{g - Qj)S(z - Zj)) /N(N - 1) (the average is on 

the ground state). Figure ||(a) shows a contour plot of P(g, z; go, zo) at various values of d, 
one per column. An electron is fixed (black bullet) at the position [go, zq) in one dot, at the 
maximum of the charge density: the contour plots of the top (bottom) row, with z — zq — 
dot 1 (z = dot 2, zo = dot 1) fixed, represent the conditional probability of measuring an 
electron in the xy plane of dot 1 (2), given a first electron is fixed in dot 1. We find that the 
three phases we discussed above clearly show characteristic structural configurations (see also 
insets in Fig. 1): (I) At small d (left column), the whole system is coherent. The electrons, 
delocalized over the dots, arrange at the vertices and the center of a regular pentagon. (II) 
At a critical value of d — 4.6 nm (center) an abrupt transition takes place. In the new 
arrangement the electrons are at the vertices of a regular hexagon. Contrary to phase I, the 
peaks in the upper and lower dots have different heights. (Ill) When d is further increased, 
the structure continuously evolves into two isolated dots (right), coupled only via Coulomb 
interaction. Three electrons in each dot sit at the vertices of two equilateral triangles rotated 
by 60 degrees. 

Next we compare our system with its classical counterpart |l3|] . To this aim, we calcu- 
lated the (circularly symmetric) radial pair correlation function g(g; z, zo) = fdgo P(go + 
g, z; go, Zq) vs g, which gives the probability of finding an electron at a distance g from an- 
other one fixed on the same or opposite dot. For a classical system at zero temperature, 
g consists of a set of (5-function peaks. Figure ||(b) shows the calculated g(g; z, Zq) for the 
same three configurations of ||(a) , together with the histograms representing the equilibrium 
configuration of the corresponding classical system |l3|,[2i|. When we fit g to a set of gaus- 
sians, also shown in Fig. ^(b), we see that the agreement between "classical" and "quantum" 
cases is remarkable in phase III and, to a minor extent, in phase I, both from the point of 
view of absolute positions and relative intensities of the peaks. In contrast, in phase II the 
peaks of the quantum system resemble the classical ones just in position, while the intensities 
are definitely different; calculations show, moreover, that the relative intensities are strongly 



4 



EUROPHYSICS LETTERS 



d-dependent. Classically, only phases I and III are ground-state configurations JL3], while the 
hexagonal structure II is a metastable state ]To| . Quantum fluctuations thus force a novel 
phase to appear. 

We interpret the sequence I — » II — » III as a 2D — > 3D — > 2D transition. In phase 
III the motion is quasi-2D, with electrons occupying degenerate S and AS states to form a 
staggered configuration: the equilibrium structure corresponds to the classical one, with inter- 
dot tunneling absent. In phase II, instead, the motion acquires an effective z-component: as d 
is decreased, electrons keep their electrostatic repulsion energy low at the expense of occupying 
AS states separated from S orbitals by an increasing energy gap. There is no classical analogue 
to this phase. When the gap becomes too large, the electrons abruptly arrange into I, all frozen 
in S orbitals. The motion now is that of a coherent quasi-2D QD, comparable to a classical 
single QD. To sum up, only in phase II electrons move in a truly 3D system: as shown in 
Fig. ||(a), the angular modulation of P is much weaker than in I and III, suggesting that the 
crystallization regime has not been reached yet (2[ |) . 

There is another, more subtle, discrepancy between quantum and classical cases, due 
to the interplay of tunneling and particle indistinguishability. In phase I classical particles 
belong either to one dot or to the other, three each for N — 6, and they can arrange in very 
asymmetric configurations in each dot, like in Fig. 1 of Ref. |L3]], provided that overall they 
form the centered pentagon which is favored by Coulomb interaction, the only energy scale 
in this case. In the quantum case, on the contrary, particles are equally delocalized on both 
dots. It simply makes no sense in this phase to assign electrons to a specific dot, since at small 
d they can easily penetrate the inter-dot barrier. Therefore, tunneling changes the physics of 
the artificial molecule and actually drives phase I. Indeed, if we suppress tunneling (e.g., via 
a sufficiently high barrier), phase III is the only one present at any value of d, while phase I 
never shows up, contrary to the classical prediction [jl3|l . 

We analyze the spectrum of neutral collective elementary excitations computing the dy- 
namical form factor S(L,oj) = |(n|p^|0}| 2 S(uj — u> n + u>o), with \n) (|0)) many-body n-th 
excited (ground) state with energy u> n (oj ), u> excitation energy, p\ = J2 m J2 g , g > c L+ m , g c ™,g' 
L-th angular component of density fluctuation. S(L, lo) assigns a weight to a charge density 
wave of angular momentum L and energy u in the excitation spectrum. Figure [|(a) shows 
S(L, lu) at low energies in the d — w space; the intensity is proportional to the radius of circles. 
Only few branches are dominant in each phase, with a characteristic value of \L\: 5 in I, 
6 or 3 in II and III, respectively. The occurrence of very few excitations with large values 
of S(L, ui) p4| is remarkable: since the ground state is formed by a linear combination of 
thousands of Slater determinants with non-negligible weight, high values of the form factor 
are due to a constructive interference effect. Indeed, the geometry of the Wigner molecule 
selects the allowed -magic- values of M and parity, and thus L = — M\ corresponding to 
the coupling between two magic states Mi, M^. Consider, e.g., the wavefunction of phase 
II with electrons at the vertices of a regular hexagon. A cyclic coordinate permutation is a 
7r/3 rotation such that ^ — > exp (7riM/3)^. In this case ^ changes sign [ p5[ , and therefore 
M = 6p + 3 (p integer) and odd parity. Our findings show that an excited magic state can be 
created adding L quanta of angular momentum to the ground state via the density fluctuation 
operator p\. The allowed values of L characterize different phases. The diamond shape of 
different |L|-branches in Fig. ^(a) originates from level crossing between magic states: uj — > 
at every M-transition in the ground state. 

Finally, we compute the "electronic" part of Raman scattering cross section |26). This 
is strictly related to S(L, oj), and is proportional to ^ n |M„o| 2 6 (lj — uj u + u>o), with M n o = 
J2 a a { a \ eiq r \P) ( n \ c a c 3\Q)> an d Q wave vector transferred in the inelastic photon scattering 
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[ 27 1 . The correlated excitation spectrum of single QDs has been already experimentally probed 
[28 1, in the regime ql « 1. Thus, having fixed a suitable value of q in a backscattering 
geometry, in Figs ||(b) and (c) we plot the cross sections at different angles 9 between q and 
the xy plane (9 — 0° and 45°, respectively). The spectra resemble the form factor S(L,ui), 
but intensities are differently modulated. For in-plane scattering [Fig. ||(b), 9 — 0°], only the 
branches with \L\ = 5 characterizing phase I survive, while all other signals are suppressed. 
The background small dots appearing in Fig. |3|(b-c) are smaller than 9 % of the maximum 
value. As 9 increases [Fig. ||(c), 9 = 45°], the branches of phase I continuously loose intensity 
as those of II and III acquire weight, until the latter show a very strong signal at 9 = 45°, 
while the former are suppressed. In I, indeed, all electrons occupy S orbitals, hence the z- 
component J dzx*(z) e lq * z Xj(z) of matrix element (ale 19 ' 7 "!/?) diminishes as long as 9 (and q z ) 
increases; the opposite holds true for II and III, since S and AS orbitals can mix. This means 
that momentum can be transferred in the z-direction only in II and III, contrary to I, where 
the system is a unique quasi-2D dot. Signals with |L| = 6 are almost invisible, but still phase 
II can be resolved from III in virtue of the characteristic energy and slope of visible branches. 
This comes from the way S and AS orbitals are differently filled in II and III, respectively, 
and hence from the different contribution of tunneling to the total energy. When 9 is further 
increased, the selection rule L — becomes effective, and eventually at 6 = 90° low-energy 
signals are suppressed. 

To summarize, we have predicted a 2D-3D-2D transition of interacting electrons in a 
double QD, accompanied by the appearance of a novel liquid-like quantum phase, in addition 
to classical configurations. The ^-dependent modulation of the Raman spectrum should allow 
to experimentally discriminate between different phases. 
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inter-dot distance d (nm) 



Fig. 1 - Ground-state energy (left axis) and in-plane average radius {q} (right axis) vs d for N = 6 
at B — 25 T. For large d- values M = 36 (not shown), twice than in a single QD with N = 3. Insets 
show the electron arrangements in the different phases. 
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in-plane radial distance p (Z units) 



Fig. 2 - (a) Contour plot of the pair correlation function P(q, z; go, zo) at various d-values. An 
electron is fixed at (qo,zq) (black bullet) corresponding to a maximum of charge density, (b) Radial 
pair correlation function g(g; z, Zo) vs g. g has been renormalized in such a way that J dgg(g; z, zq) = 
1. Lengths are in units of £ = 7.20 nm [£ = (Ti/m*^l) 1 ^ 2 characteristic in-plane radius]. Insets 
schematically show the spatial arrangement of electrons projected on the xy plane. 
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inter-dot distance d (nm) 



Fig. 3 - (a) Dynamical form factor S(L, u>) in the d — u> plane, for different angular momentum 
components L, indicated in parentheses. The radius of each point is proportional to its intensity, (b) 
Analogous plot for the Raman scattering cross section, at an angle 6 — 0° between the xy plane and 
the transferred momentum q (q = 2 ■ 10 6 cm -1 ), (c) Same as (b), with = 45°. 



